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Abstract. We consider quantum phase transitions out of topological Mott insulators 

in which the ground state of the fractionalized excitations (fermionic spinons) is 

topologically non-trivial. The spinons in topological Mott insulators are coupled to 

i-rt ' an emergent compact U(l) gauge field with a so-called "axion" term. We study 

fjj . the confinement transitions from the topological Mott insulator to broken symmetry 

phases, which may occur via the condensation of dyons. Dyons carry both "electric" 
and "magnetic" charges, and arise naturally in this system because the monopolcs 
of the emergent U(l) gauge theory acquires gauge charge due to the axion term. 
It is shown that the dyon condensate, in general, induces simultaneous current and 

£f) ■ bond orders. To demonstrate this, we study the confined phase of the topological 

Mott insulator on the cubic lattice. When the magnetic transition is driven by dyon 
condensation, we identify the bond order as valence bond solid order and the current 
order as scalar spin chirality order. Hence, the confined phase of the topological Mott 
insulator is an exotic phase where the scalar spin chirality and the valence bond order 

£S| ■ coexist and appear via a single transition. We discuss implications of our results for 

generic models of topological Mott insulators. 
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1. Introduction 

One of the most notable progresses in modern condensed matter physics is the 
development of the concept of "topological order". Topologically ordered phases 
cannot be understood in terms of order parameters, which are the bases of Landau- 
Ginzburg theories of broken-symmetry phases. Standard examples of the topological 
order include quantum Hall phases [TJ |2] and time-reversal symmetric topological 
insulators [31 IU EJ EJ [7] (we will abbreviate 'time- reversal symmetric topological 
insulator' as 'topological insulator'). Topological insulators are insulating in bulk but 
have gapless edge or surface states which are 'immune' to the opening up of a gap as long 
as the time-reversal symmetry is respected. Though the protected gapless edge/surface 
state itself is quite interesting, it is notable that the topological band insulator can 
provide a condensed-matter example of axionic electromagnet ism P, EJ [TO] 

U = ^- 2 e^ Xp F, u F Xp , (1) 

when the surface is gapped by the time-reversal symmetry breaking perturbations. In 
condensed matter physics langauge, this axionic term can be understood as a half- 
integral surface quantum Hall effect [9] . This quantized anomalous Hall effect and other 
interesting surface and bulk physics of the topological band insulator have attracted 
both the theoretical and experimental interest [HI H21 [131 [IS HSl HSl 113 HBl HSl I2DI 121] ■ 

Axionic electromagnetism, especially in the presence of a dynamical gauge field, 
has a long history in high-energy physics. The interest in the axion term stems from the 
CP- violation problem [22| |23| 121] and confinement of quantum chromodynamics [25J . 
Though the axion term equation [1] is a total derivative, 't Hooft [25] showed that the 
axion term is capable of changing the nature of the instantons and the confined phases 
in non-abelian gauge theories. The phase where the dyon is condensed in the presence of 
the axion term has a qualitatively different structure from the ordinary confined phase; 
this difference is traced back to the fact that the monopole receives the gauge charge due 
to the axion term (this effect is sometimes called "Witten effect" [26J). Similar behavior 
can emerge in abelian gauge theories on the lattice; Cardy [27J [28] studied the phase 
diagram of the gauge theory with the axion term and demonstrated oblique confinement 
in the discrete Z p gauge theory and the compact U(l) gauge theory. The structure of 
the phase diagram and the related SL(2; Z) symmetry are well studied [271 EH1 129] . 

As we now have a condensed matter system exhibiting axion electromagnetism, 
we are led to the natural question: could we find an example where the axion term 
changes the nature of the phases and phase transitions in the condensed matter system? 
The monopole in the presence of the axion term becomes a dyon with the fractional 
charge [251 ESI EQ] , and the structure of the confined phase should reflect this fractional 
charge of the dyon. In this paper, we will seek to find a condensed matter example of this 
novel confined phase which reflects the axionic electromagnetism and dyon condensation. 
To study the confinement, it is important to allow magnetic monopoles in the excitation 
spectrum. At the same time, the compactness of the U(l) gauge theory is important 
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for the existence of monopoles. Hence, to study the confined phase with Witten effect, 
we need both a compact U(l) gauge theory (or non-Abelian gauge theory) and the 
axion term. Clearly, the usual electronic topological band insulators with the non- 
compact "external" electromagnetic fields cannot provide such an example. In the 
condensed matter system, the compact U{\) gauge theory on the lattice often emerges 
as a consequence of "fractionalization" of the electron. It has been suggested that 
the so-called topological Mott insulator, where the fractionalized excitations or spinons 
possess topological band structure, satisfies both requirements, i.e., fractionalization 
and topological band structure [31] . 

There may be two different routes to obtain topological Mott insulators. First, 
one could start from an electronic topological band insulator [31]. Upon increasing the 
on-site Hubbard interaction, there may be a transition to a spin liquid Mott insulator. 
Such a transition can conveniently be described by the framework of the slave-rotor 
field theory, where the electron operator is written as the product of a fermonic spinon 
operator and a bosonic rotor operator with proper constraint on the Hilbert space. Thus 
the low energy effective field theory naturally becomes a U(l) gauge theory, where the 
fermonic spinons and rotor bosons are coupled to a U(l) gauge field. Here the fermonic 
spinons carry the spin- 1/2 quantum number and the rotor bosons the electric charge 
of the electron. The transition from the spin liquid to the electronic topological band 
insulator can be regarded as the bose condensation transition of the rotor bosons. In the 
spin liquid phase, the charge-carrying rotor bosons are gapped and the fermonic spinons 
would inherit the topological band structure of the electrons. These fermonic spinons 
are in turn coupled to an emergent U(l) gauge field as described above. When the 
fermionic spinons are integrated out, the non-trivial topology of the spinon spectra leads 
to an axion term for the emergent U(l) gauge field. In the second case, one may start 
from spin models [32J with frustrated antiferromagnetic and ferromagnetic exchange 
couplings. It was recently shown that slave-fermion theory of such spin models can 
support a U(l) spin liquid, where the fermonic spinons acquire an emergent spin-orbit 
coupling. In some frustrated lattices, this leads to the topological band structure of the 
spinons in analogy to the non-interacting electronic topological band insulators. In this 
case, the non-trivial topology of the spinon spectra is emergent and not inherited from 
the electrons. Similarly to the first case, the axion term would arise upon integrating 
out the spinons in the bulk (we refer the readers to the reference [311 [32] for the detailed 
discussion of the transition and phases) 

In order to gain some insight as to the dyon condensate in the topological Mott 
insulator, we start by briefly reviewing the confined phase and the Higgs phase [33j [34] , 
driven from a "bosonic" U(l) spin liquid by condensing the monopole and the "bosonic" 
spinon. The condensation of bosonic spinons leads to the conventional Neel phase. 
On the other hand, we obtain the Valence Bond Solid (VBS) state if we condense 
the monopoles in the U(l) spin liquid in (3 + 1)D, which can be understood as the 
confinement of electric fluxes. In the VBS state, we have confined spinons and no 
"photon". As there are significant connections between the ordinary confined phase 
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and the oblique confined phase, we discuss first the nature of the ordinary confined 
phases. First, the non-trivial spatial patterns of VBS originate from the crystal 
momentum carried by the monopole operators. The crystal momentum can be traced 
back to the nontrivial Berry phase due to the background staggered U(l) gauge charge. 
The staggered U(l) gauge charges induce the background "electric" fluxes connecting 
charges. When a monopole hops around the dual lattice, it sees the "electric" fluxes and 
feels the Aharanov-Bohm phase. This is exactly dual to the standard Aharanov-Bohm 
phase for electric charge circulating around magnetic flux. Due to the Berry phase for 
the monopole, the monopole tends to condense at a finite momentum. Condensation of a 
monopole at a finite momentum implies that there is a non-trivial spatial pattern of the 
monopole current [331 US]- Then, the monopole current will induce the "electric" fluxes 
just like an electric current inducing a magnetic field. As the valence bond operator is 
proportional to the "electric" flux, we end up with patterns of valence bonds. 

Some of the above discussions for the standard confined phase can be applied to the 
confined phase driven by dyon condensation in the topological Mott insulators. First 
of all, the dyon will experience a Berry phase if there is a nontrivial background charge 
pattern because the dyon carries both the "magnetic" and "electric" charges. It will be 
shown in the main text that there exists another contribution to the Berry phase due to 
the axion term in the topological Mott insulator. Secondly, when the dyon with a non- 
zero crystal momentum condenses, the system supports both the "electric" current and 
"magnetic" monopole current. The dyon current then will induce the 'magnetic" fluxes 
as well as the "electric" flux. The "magnetic" flux should be naturally accompanied by 
the gauge-charge current in the oblique confined phase. Hence, we obtain a nontrivial 
phase where the gauge-charge current order coexists with the bond order when the dyon 
condenses. 

There are a number of subtle issues in dyon condensation and the physical 
interpretation of such phenomena. First, it is known that the statistics of the dyons 
may not be trivial [351 [36l [3TI [381 EHl SQl SU S21 H3] . We show that there is no 
statistical transmutation for the dyons in the topological Mott insulator. Thus the 
dyons are bosons as far as the monopoles are bosonic [381 H2] ■ Secondly, we argue that 
the monopoles in topological Mott insulators have as their quantum numbers only the 
U(l) gauge charge and crystal momenta. (From the discussion of the usual topological 
band insulator P [10], we physically understand the gauge charge of a monopole as 
a 'polarization' charge so that the dyon is not a bound 'particle' to the monopole. 
This understanding of the charge of the dyon perhaps implies that there will not be a 
statistical transmutation as the charge is more like a polarization 'cloud', not a particle). 
The actual pattern of the dyon condensate depend on the details of the microscopic 
physics, i.e., the lattice structure and the background charges. In this work, we discuss 
the dyon condensation patterns and their physical implications on the cubic lattice. 

The rest of the paper is organized as follows. In section II, the field theory of 
dyons on the dual space/lattice is derived starting from the U(l) gauge theory with 
axion term on the direct space/lattice. We discuss possible dyon condensation patterns 
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and resulting confined phases in section III. We summarize our results and conclude in 
section IV. 

2. The field theory of dyons on the dual lattice 

We consider fermionic spinons on a cubic lattice for simplicity to develop the basic idea 
of dyon condensation. In the following section we connect the general considerations in 
this section to possible phases proximate to the topological Mott insulator state. 

In order to derive the action for the dyons (or monopoles), we first need to integrate 
out the gapped matter fields in the system and obtain an effective gauge theory. As the 
fermionic spinons are in the topological insulator state, the topological nature of the 
spinon states leads to an axion term in the effective U(l) gauge theory as follows: 

La = Lq ED + Lq (2) 

The first term Lqed describes lattice quantum electrodynamics in (3 + 1)D in the 
presence of the staggered gauge charge 

Lqed = ^~ 2 {d,A u - d v A,f + i£ A T (r)N(-iy + « +z , (3) 

where N is the magnitude of the staggered gauge charge for the U(l) gauge field 
residing on the two sublattices of the (bipartite) cubic lattice. That is, there exist 
background gauge charges (— \) x+y+z N at the lattice points (x,y,z). As will become 
clear later, these background charges are the source of the Berry phases experienced 
by monopoles (or dyons). The existence of the background charges on the lattice 
for the monopoles/dyons can be seen from the following argument. If there were 
no background charge, the dispersion of monopoles/dyons on the dual lattice would 
have the minimum at the zero crystal momentum. The confinement transition via the 
condensation of monopoles/dyons would occur at the zero crystal momentum. This 
would suggest that the resulting paramagnetic insulator (with a charge gap) would not 
break any translational symmetry and at the same time there would be no ground state 
degeneracy. When the number of spins/fermions in the unit cell is odd, this is not 
possible because the translationally invariant paramagnetic ground states in such cases 
can only have either a gapped phase with a ground state degeneracy or a gapless phase 
via the Hasting's theorem jH] . Thus it is clear that there should be a background charge 
for the monopoles/dyons in the case of the odd-number of spins per unit cell. When the 
number of spins/fermions per unit cell is even, which is the case in the topological Mott 
insulator, this is not guaranteed as the Hastings' theorem cannot directly be applied. 
Nevertheless, in light of the discussion for the case of the odd-number of spin per unit 
cell, it is conceivable that there should at least exist the situations where there is a 
non-trivial background charge for the monopoles/dyons. This is the case that we study 
in our work. 
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The second term Lq of equation [2] is the axion term for the emergent compact U{1) 
gauge field, arising from the topological properties of the spinons: 

° e^d p A v d x A p . (4) 



3tt 2 



S = / drd 3 r — m ^—- m + i(Xg + — a°)m" . (6) 



To perform a duality transformation on the corresponding action, 

S A = J drd?v^- 2 {d,A v - d v A p f + ±^e^d,A v d x A p + i ]T A T (r)N(-iy +y+z , (5) 

we introduce the monopole current m p = j-e^ vXp d v S\p, where S\ p G 2ix7h so that 
m p e Z. Introduction of the monopole field m p modifies Sa such that d p A v — d p A u — > 
d p A v — d p A u + S pv . Since the action is quadratic in A p , we can formally integrate out 
A p . In order to do so, we first solve the Gauss law constraint for the "electric" field of 
the U(l) gauge theory, namely V ■ E = N(—l) x+y+z , by breaking the gauge field A p into 
the fluctuating part a p (which is free of the constraint) and the static field a that is 
responsible for the staggered gauge charge. That is, we write A p = a p + a° p and integrate 
out a p . This generates only two terms that are simple and intuitively understandable: 

The first term describes the Coulomb interaction ~ —^between dyons with the modified 
strength G 2 = (f^) 2 + %- due to the "electric" charge g/2 = g x ^ with 9 = n and 
"magnetic" charge 2ir. The second term represents the Berry phase experienced by the 
dyons due to the staggered QED charge V • E — N(—l) x+y+z . Here, Xq is the dual 
£7(1) gauge field for a . This term has the form of the minimal couplings of m p to a 
and X®, and simply dictates that the monopole carries the U(l) gauge charge 9/2n and 
the unit magnetic charge 2it as expected from the axion term. 

We now derive the dyon action in equation [6] rigorously, and also show that there 
is no long-range statistical interaction between the dyons, that is, that the statistics 
of the monopoles does not experience any statistical transmutation when they are 
converted to dyons. We begin with the Lagrangian of the action Sa in equation [5] 
with A p = a p + a° p where a° p represents the static configuration due to the background 
gauge charge N(—l) x+y+z . Splitting the Lagrangian in terms of a and a p , and using 
fp, v = d p a u — d v a p and /° = dpO? v — d u a° p , the resulting Lagrangian can be written as 

L = Li[a^, S pu ] + L 2 [a° p , S pu ] + L 3 [a p , oj]. (7) 

The Lagrangian for the purely fluctuating part, L 1; consists of the standard QED 4 
supplemented by the axion term in the presence of the monopole currents. 

1 Q 

Li[a p , S pv \ = ——^{jfj,u + s pu ) + 2 £ Jnub\p+ 
9 9 

c u pu u \p ~T „„ 2 C JpuJXp 



32tt 2 "" Ap 32vr 2 

The other parts of the Lagrangian are given by 

L *K s ^\ = ^O 2 + « + ^C^ + tJ^ a %°a, (9) 
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where rf" = N(—l) x+y+z 5^' T encoding the staggered QED charge, and the mixed 
contributions 

LsK, oj] = ^ 2 fU^ + 4^ £ " VXPf ^ hp - (10) 

Up to this point, there have been no approximations. We now transform each term into 
a desired form via approximations justified in the limit of the large mass gaps for the 
monopole and the spinon. We begin by dropping L^a^a ^] entirely using dpf^ v ~ 
and e tJlUXp d u fx p = 0. The first condition d^fp U ~ is satisfied if the gap for "charged" 
excitations is large enough that there is no free current giving rise to a field strength. 
In the Coulomb phase, every matter field (including monopoles and dyons) is gapped 
so this approximation is justified. For L 2 , we solve ^r(/^) 2 + ir/^a ^ first to fix a° and 
end up with the following Lagrangian 

i n 

Zg 107T qjn 

where the equality /° = "^e^xpdxX (X° ~ dual U(l) gauge field) and m M = 
j^^puXpdySxp are used. Thus, L 2 contains all the important Berry phase terms for 
the monopole operator. Now we integrate out a M in L\ to investigate the Coulomb 
interaction between monopoles and the statistical interaction. We choose the radiation 
gauge for the fluctuating U(l) gauge field a M (i.e., we use the "photon" propagator y ^2) 
to obtain 

= 2 2-7T m ' il Tf m ' ii ~^~ A^Z ^ A"' "*" 2VfJ,b al j,TJ2-O u bau) 

a a 1 

+ 77?; O^ bnyDXp + - ^OyDny-^-TE O a OR p , {*-2) 

32tt z 8n z o z 

Straightforward calculation reveals that the first line in the above L\ becomes the 
Coulomb interaction [28j EH] between dyons and the second line exactly vanishes, i.e., 
-^e' MuXp S^ u Sxp + -^d v S ^jpE^PdaS p p = up to the boundary terms. We are discussing 
the bulk transition without boundary, hence this term can be dropped out completely. 
In fact, this term corresponds to the ^-dependent statistical interaction between the 
dyons as discussed in A. Goldhaber et.al. [3H] Hence, the dyons do not receive additional 
statistical interaction among themselves. This suggests that the dyon in the topological 
Mott insulator is bosonic and can be condensed. Using 6 = 71 and collecting the resulting 
terms in L\ + L 2 + L 3 , we arrive at the effective action for the dyons 

S = J olrdh f ?! m M_L m/l + i(X£ + \al)m»\ , (13) 

as advertised before. To proceed further, we introduce the soft-spin operator ijy ~ 
exp(— i<f>) for the dyon creation and put the lattice structure back to obtain the 
Hamiltonian on the dual lattice. 

H = -tJ2 4 exp(-t(X* R ' + ^a° RRI + GL RRI ))^ RI + h.c, (14) 

R,R' 
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where we introduce an auxiliary fluctuating U(l) gauge field L RR/ to encode the Coulomb 
interaction with the coupling constant G = a/ g 2 /4 + 47r 2 /g 2 , and R denotes a dual 
lattice site. 

Some comments on this Hamiltonian are in order. First of all, we need to be 
careful in defining a RR , on the dual link as a was originally meaningful on the links of 
the direct lattice. To resolve this issue, we need to introduce a short-ranged function 
Fffi which maps the direct link to the dual link. This defines a RR , = J2<rr>> ^rr^ a rr'- 
The appearance of this short-ranged function Fffi can be traced back to the form 
of the axion term which connects f^ u to f\ p by s^ vXp and we are forced to introduce 
such a function f27\ |2"B"] if we wish to study the axion term on the lattice. Secondly, 
we have introduced an auxiliary non-compact U(l) gauge field L RR > which differs from 
a RR i and X RR >. In fact, this fluctuating gauge field corresponds to the U(l) gauge field 
which is rotated in the electric-magnetic charge. This gauge field is introduced only 
for convenience as there is no corresponding physical charge. The two physical charges 
in the dyon theory are the "electric" charge and "magnetic" charge of the compact 
£7(1) gauge theory A^. However, dyons carry both charges, and the effective Coulomb 
interaction between the dyons is the sum of the Coulomb interactions of "electric" 
charges and "magnetic" charges. Because of this, the dyons are seen to interact with 
the gauge charge G = a/ g 2 /4 + Arc 2 / ' g 2 . This interaction can be written compactly 
at the cost of introduction of the auxiliary U(l) gauge field L RR /. As the low-energy 
physics is determinded by the static component, we can ignore this to study the mean- 
field physics [31] . In the case when a° rr , = can be choosen (no flux threading the direct 
plaquette), we need to solve the following single-particle hopping problem 

H = -tj^ ^exp(-zX i?i? ')^' + h.c. (15) 

R,R' 

This Hamiltonian leads to a set of degenerate energy minima for ip R . This degeneracy 
will be lifted by the non-linear terms mixing these degenerate minima. This is the same 
procedure as used in studies of monopole condensation and the resulting confinement 
problem. 

There is, however, an important difference from the monopole condensate problem: 
the number of dyon species. We now argue that it is enough to consider one doublet 
of dyon operators in the low energy limit. To see this, we first note that the monopole 
of the strength 2ir can collect the U(l) charge n + 1/2, n e Z when the axion angle 
9 = it. Then each dyon will have the self-energy E ~ (n + l/2) 2 g 2 + 4ir 2 /g 2 . Thus, the 
dyon of the monopole strength 2n has two degenerate states with the charge 1/2 and 
— 1/2, representing the lowest energy dyons. We label ipi and ip2 as the dyons of the 
magnetic and electric charge content (2ir, 1/2) and (27r, —1/2). These two fields (ipi,ip2) 
are connected to each other by the time reversal symmetry as the time-reversal operation 
T will flip only the strength of the magnetic field, i.e., T[ipi] = ip 2 an d T[tp 2 ] — iA where 
f operation (particle to anti-particle operation) flips both of the magnetic and electric 
charge. 
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The other way to understand the nature of the doublet is to consider the time- 
reversal symmetry of the spectrum of the dyon fields. We have the time-reversal 
symmetry before we condense the dyons. As we have a dyon ipi of (2ir, 1/2), we need 
to have a time-reversal partner of this state, i.e., (— 2n, 1/2) which is nothing but the 
state ipl- Hence, we see that when ip% (or ip 2 ) is condensed, we enter a phase where the 
time-reversal symmetry is broken. We note that even though the above derivation is for 
the U{1) gauge theory, similar calculation for Z p lattice gauge theory should be possible 
and indeed has been done by J. Cardy [271 EH] without the background charges. 

3. Dyon condensation and resulting broken symmetry phases 

In this section, we study the broken symmetry phases that arise when the condensation 
of dyons occurs in the topological Mott insulator. As the precise ordering depends on the 
lattice structure, we specifically work on the cubic lattice. To determine the patterns 
of the symmetry breaking explicitly, we work out the lattice symmetry of the dyon 
fields (more precisely, we determine the projective symmetry group [451 HEJ H7] of dyon 
operators) and solve the dyon hopping problem. The effect of lattice symmetries (at 
least translations and rotations) on the monopole operator has been thoroughly studied 
in the literature and we build on those results [34J. We consider the dyon hopping 
problem on the dual lattice defined as 

H = -t J2 ^a, R eM-iXo R ')^a,R> + h.c, (16) 

a=l,2;R,R' 

where we choose Xq R that is consistent with the background charge (— \y x+ y +z . 
Although here we work on the cubic lattice with the staggered background charge N = 1 
in equation [THl it should be straightforward to generalize this for different bipartite 
lattices and larger N G Z. For example, a larger staggered charge N = 2 on the cubic 
lattice has been studied for the ordinary confined phase, and N = 2 background charges 
were found to induce similar patterns as the case of N = 1 [33] . 

In the case of AT = 1 [34], we have two minima for each ip a where the dyon 
operator takes its minimum kinetic energy. We label these two minima by o = 1,2 
and hence we have a total of four operators ip a ,a, a = 1,2, a = 1,2. The corresponding 
eigenfunctions [34] have the forms: 



1 + (V3-V2> 



ITTZ 



Vv 



2(3 - VQ) 
(V3-V2)e 



te 



IT if 



(17) 



; 2(3 - v^) ^ 

Hence, we introduce the dyon operators as follows. 
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The operations of the spatial symmetries on ^/ aa are worked out in the paper (31] , and 
quoting their results, we get 

r«:(*a,i.*^)->(*;i,-*y 

T M :{*a,i,*afl)->(Kfl,K,i) (19) 

Ry : (*«,!, ^M "► ( ^ °' 2 , *^= Q ' 2 )- 

Note that the spatial symmetries do not mix the dyons with different U(l) gauge charges, 
i.e., the translations and rotations do not change the index a = 1,2. To change a, 
we need to consider the time-reversal operation on \l/ ajCr because T[^ lo -] = ^2 a an d 
r[ 1 i r 2,(r] — ^lcr- The time-reversal operation does not change the index a because (1) 
the dyon operator in the dual theory does not change the sign (i.e., T : tRRi — > tRRi), and 
hence no change is involved for the minima of the kinetic energy of the dyon operators 
(2) the minima of the kinetic energy of the dyon operators are at the time-reversal 
symmetric momentum. Hence, the spatial symmetry is completely determined by the 
actions on a. We are now ready to introduce the following two vector order parameters 

E r = 

B r = , (2U) 

for r = x,y,z. For V = E,B, we have the following transformation rules 

T* : (V x ,V v ,V t ) ^ (-V x ,V y ,V z ) 

T y :(V x ,V y ,V z )^(V x ,-V y ,V z ) 

T z :(V x ,V y} V z )^(V x ,V y: -V z ) 

R z :(V x ,V v ,V z )-^(V y ,V x ,V z ) 

R y :(V x ,V y ,V z )^(V z ,V y ,V x ). (21) 

The vectors V = E,B follow the same spatial symmetry transformation rules because 
the spatial symmetry is fully determined by the actions on a in \1/ Q)CT , a = 1, 2, a = 1, 2. 
In addition to the lattice symmetry, we also need to introduce the time reversal 
symmetry operations: 

T:(E,B) = (E,-B). (22) 

In the end, we have two vector operators E and B. What phases do these two vector 
operators E and B represent? It can be shown that the field E behaves like staggered 
bond order (equivalently, electric flux) and B behaves like staggered current order 
(magnetic flux). 

If the confinement induces ordering of the magnetic degrees of freedom, then the 
dual operators E and B should have representations in terms of the spin operators 
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S [32] • In the previous studies, it is already noticed that E represents the VBS-like 
order parameters [331 EI]- More explicitly, E x ~ (— l) x < S r ■ S r+X > and there are 
analogous identifications for E y and E z . This VBS corresponds to "confined electric 
flux". 

It is not difficult to see that the magnetic flux expected in the dyon-condensed 
phase is the field B. This magnetic flux should be sourced by the spinon current living 
on the plane perpendicular to the B. In turn, the spinon current can be related to 
the scalar spin chirality. These facts, hand-in-hand, lead us to the conclusion that the 
scalar spin chirality can be connected to the vector field B. This intuition is supported 
by explicit construction of appropriate combinations of the scalar spin chirality, which 
follow the same transformation rules as B. 

B x = (-l) x M yz , B y = (-l) y M zx , B z = (-l) z M xy (23) 

where M y - is the scalar spin chirality defined on the r/'-plane. For example, M xy involves 
four spins of the plaquette on the direct lattice. We label four sites 1 = (x,y,z), 
2 = (x + l,y,z), 3 = (x + 1, y + 1, z) and 4 = (x, y + 1, z) in the x?/-plane. Then, M xy 
is defined uniformly on every plaquette of the xy-plane as 

M xy =<S 1 -S 2 x S 3 > + <S 2 -S 3 x S 4 > + <S 3 -S 4 x S 1 > + <S 4 -S 1 xS 2 > . (24) 

It can be shown that (— l) z M xy and B z follow the same transformation rules under the 
lattice symmetries and time-reversal symmetry operations. Analogous expressions hold 
for B y and B z (as noted before, the scalar spin chirality lies in the plane perpendicular 
to the vector B). Thus, we confirm that the scalar spin chirality is induced in the dyon 
condensed phase. Now we discuss where the scalar spin chirality order comes from. The 
current order B originates from the 'mixing' of the two dyon species with different U(l) 
gauge charges. This order parameter would have not been present if there were only 
one monopole species as in the usual confinement problem. In the case of the dyon 
condensate, both of the electric flux and magnetic flux should be confined at the same 
time. This implies that there should be sources of the magnetic flux as well as the 
electric flux. We can easily anticipate that the magnetic flux (~ B) can arise from the 
spinon current (~ M y ) via the elementary Biot-Savart Law. 

With these identifications between the dyon operators and the corresponding 
order parameters, we can now ask what kind of phases will eventually appear. This 
requires the studies of the quartic order and higher order terms in the Landau- Ginzburg 
theory [33l [33] written in terms of the dyon fields \1/^ = (^ a ,i, ^0,2), a = 1,2. We will not 
attempt to solve the full phase diagram but will point out a few representative confined 
phases from the symmetry-allowed forms of the Landau-Ginzburg theory 

L = Luni^x) + L kin ($ 2 ) + L int ($i, \J/ 2 ). (25) 

Here, L kin represents the kinetic terms in equation [15j We expand L int = L 2 + L 4 + ■ ■ ■ 
with L n ~ 0(^ n ),n G Z. By requiring the lattice symmetry and the time-reversal 
symmetry, we identify L 2 = r|\I/i| 2 + r|\E f 2| 2 and 

L 4 = u (#} ■ ^i) 2 + u(#| • ^ 2 ) 2 + vl^l 2 !^! 2 + w\$! ■ $ 2 \ 2 . (26) 
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Up to this quartic order, we notice that there are several interesting possibilities. First, 
< ^x >= and < ^2 >= for r > 0: This phase is the topological Mott insulator 
phase. The order parameters are E = and B = 0. Secondly, when the sign of the 
mass term for \l/ a is negative r < in L 2 , two dyon fields \I/i and \I/ 2 tend to condense. 
We assume u > in L4. There are four following possibilities of L4. 

(1) If w < and 2u > v > —2m, we have < ^1 >= pexp(i(j)i),< \J/ 2 >= P*exp(z0 2 ) 
where p is a real two-component spinor. Then, the order parameter E oc (fop 7^ 
but B = 0. The two dyon fields condense at the equal amplitude in a way to respect 
the time-reversal symmetry. Hence B = is inevitable, and this phase is the ordinary 
confined phase. 

(2) If w > and v > 2u, then we have < ^1 >^ but < \t> 2 >= 0. Here, E and B 
acquire non-zero expectation values (with E || B) at the same time, and in this phase, 
the scalar spin chirality and the VBS order coexist. Furthermore, the time-reversal 
symmetry is spontaneously broken. 

— * — * 

(3) If w > and 2u > v + w > —2u, we have < \l/i >= p\ exp^^), < \1>2 >= P*2 ex p(^02) 
where p\ and p<i are the real two-component spinors with p\ • p 2 = 0. In this confined 
phase, the order parameters E = but B ^ 0. Thus, this phase corresponds to the 
phase where only the scalar spin chirality orders. Here, the two dyon fields condense 
at the equal amplitude in a way to break the time-reversal symmetry maximally. This 
phase should be considered as an "opposite" of the phase (1) where only E acquires 
non-zero expectation value. 

(4) If w < and v + w > 2u, the only one of the two dyon fields condenses. The phase 
is identical to the phase (2) and there are two orders present in the phase. 

The detailed spatial patterns of the order parameters can depend on the higher- 
order terms in L int which we do not consider (in fact, the pattern depends on the 
~ 0(x 8 ) in the Landau-Ginzburg theory, see for example Motrunich and Senthil [34]). 

We would like to point out that VBS order can be detected by measuring 
the triplon spectrum in inelastic-neutron scattering |48j combined with the magnetic 
susceptibilities measurements [10]. The scalar spin-chirality may be detected via Raman 
scattering J50J [51] or resonant inelastic x-ray scattering(RIXS) [52]. 

4. Conclusion and outlook 

In this paper, we have studied the condensation of dyons and the corresponding confined 
phases in topological Mott insulators. Using the effective field theory and projective 
symmetry analysis on the dyon operators, we identified possible broken symmetry phases 
resulting from such confinement. Because the dyons are monopoles endowed with gauge 
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charge, the confined phase can induce at least two different orders: bond order and 
current order. In the case when the axion angle 9 is not n, we might need to consider 
more species of the dyons. However, the dyon carries the "electric" and "magnetic" 
charges as far as 9 ^ 0, and this implies that the confined phase in the non-zero axion 
angle 9 will be a coexisting phase of the bond order and the current order if 9 ^ 0. In this 
paper, we concentrate on possible confined phases on the cubic lattice, but our results 
should be straightforwardly generalizable to other bipartite lattices, where the precise 
spatial patterns of bond and current orders would depend on the lattice structure. 

For the physical spin/magnetic degree of freedom, the physical order parameters in 
the confined phase were shown to be the scalar spin chirality order and the VBS order. 
It should be emphasized that the coexistence of two order parameters discussed above 
is unique to the confined phase of the topological Mott insulators and does not arise in 
ordinary confinement physics in condensed matter systems (monopoles are not dressed 
with any quantum numbers in the ordinary confined phase [33], [31] ) . Notice also that the 
coexistence of the VBS and scalar spin chirality order is difficult to obtain in a Ginzburg- 
Landau theory framework in ordinary condensed matter systems as this would require 
extremely high-order coupling between spin operators while this can occur naturally if 
one starts from topological Mott insulators. Thus the emergence of the coexisting scalar 
spin chirality and VBS orders via a second order transition from a quantum paramagnet 
would be strong evidence that such a quantum paramagnet was in fact the topological 
Mott insulator. 
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